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We propose a model for motor proteins based on a hierarchical Hamiltonian that we have pre- 
viously introduced to describe protein folding. The proposed motor model has high efficiency and 
is consistent with a linear load-velocity response. The main improvement with respect to previ- 
ous models is that this description suggests a connection between folding and function of allosteric 
proteins. 

An important class of proteins are alosteric: part of the molecule undergoes large conformational changes as a result 
of a chemical reaction which takes place in another region of the protein. One example of this is motor proteins which 
OS ' through a series of conformational changes perform directed motion using the chemical energy gained ATP hydrolysis. 
Physicists have proposed that this may be viewed as a ratchet-like mechanism, similar to Feynman's famous ratchet 
JjJ but with the addition that ATP binding/hydrolysis has the ability to turn this ratchet on and off. Several models 
have been proposed along this lines, see [@-|§]- 

The energy associated with the chemical cycle of real motor proteins is comparable to the total binding energy of 
a protein. This leads us to speculate that the directed motion associated to motor action may be connected to the 
pathway for protein folding. Thus, we propose to describe the motor variables in terms of a pathway parametrization 
■ of protein folding using three variables (fo, <fi and (f2, that each takes value or 1. These variables, which 

label the conformations visited during the motor cycle, also label subsequent states along the folding pathway of the 
protein. Thus we propose that the energy function of the motor protein in the absence of ATP and external load is 
described by the Hamiltonian JToj ] : 
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The ground state of H is (ipo, ipi, (p^) = (1, 1, 1). At low temperatures T the variables will tend to freeze into the 
Q\ ' ground state, whereas they melts for T > T c = 1/ ln(2). In order to take into account ATP hydrolysis and an external 
OS ' load we suggest the following extension of the energy function: 
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where [ATP] represents the effects of ATP on the chemical potential of the molecule. [ATP] can cycle between two 
values, say and — A where A is the energy associated to hydrolysis, and F is the external load on the motor. With 
this structure we single out one variable tpo to represent ATP binding and hydrolysis while the two other variables 
O ■ represent events associated to attachment (ifi) and mechanical motion ((^2) of the motor. 

The dynamics of the motor is generated by switching ATP between —A and 0. When ATP is zero, po is driven to 
• , (p Q = 1 by the first term in the energy function. Subsequently this sets an ordering where typically first the variable 
(pi, then the variable <p2 takes the value 1. Thus first the "leg" binds to the substrate, then it performs a step (or 
"power stroke") tp2 = — ► 1, and it advances forward by one step. In case it by chance do the opposite, i.e. if tp2 
moves before ipi it represent the event when the stroke takes place before the binding, and we assign it a one step 
backwards movement. When ((^1,(^2) = (1, 1) we initiate the chemical switch [ATP] = — > —A which subsequently 
forces tpo — » 0, and accordingly release the binding constraints of (pi,<p2- The ordering of the melting of these variables 
is random. The asymmetry in the freezing — where the process typically proceed in the order ipo — ► 1, ipi — > 1 and 
tf2 — » 1 — and the melting — which can proceed in any order, we refer to as the mirror effect fic|| , and lies at the 
heart of the functioning of this motor. The total move table is 
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In terms of the normal nomenclature in the literature on motor molecules, the moves of Eq.(3) represent respectively 
the binding and the working stroke whereas the moves in Eq. (5) undo these. The moves of Eq. (6) then represent 
the optimal sequence of unbinding and recovery. 

In order to study the model numerically, we have chosen to use the following dynamics: In each update, either 
zero or one of the ip- variables changes value. The probability that a given change is made (either (<po> <£i> ^2) —> 
(1 — tpo, ipi, 1P2), (<£>0) 1 — ipi, (p%), {po, <Pi, 1 — <Pz), or (cpo, tpi, 1P2)) is proportional to the Boltzmann factor associated 
with the configuration the model is entering into, 

P( Si -» fl/ ) oc e - a ^'^ T . (7) 

where Sj and sy are initial and final states respectively. As the variables are discrete, it is not possible to assign a 
physical time to the process. However, the amount of ATP that is hydrolyzed (i.e., the number of times A takes the 
value —a) is a natural clock in this process. 

A more realistic dynamics for this model is based on constructing a Langevin equation from the Hamiltonian (|^). 
Continuous variables may then be used and a connection with "real time" may be made. The Langevin equation 
reads 

d<Pi dHi 

where fa are mobilities, %(t) thermal noise and the variables ipi now has to be continuously varying Jl^ ]. To obtain 
the real cycle times we in addition have to introduce the residence times in the different states. We note that if the 
limiting step of the dynamics is the ATP binding, i.e. at very low ATP concentration, then the efficiency curves given 
by the phase space method also represent real time velocity curves up to timescale that is given by ATP hydrolysis 
rate. In the opposite limit, the Langevin dynamics captures the timescale of the motion through the coefficients fa. 
Experiments show a roughly linear load- velocity characteristics Jl3| . We can obtain a linear dependence with this 
Langevin dynamics for certain realizations of the force but this feature appears not to be robust with respect to the 
functional form of the load. 

In Figure |I| we show the average position of the motor as function of the number of times an ATP molecule was 
hydrolyzed, i.e., that the ATP was switched from to —A. In the figure we show two cases, both for A = 4, one at a 
low temperature T = 0.25 (the upper curve) and one at a high temperature, T — 5. In each case the corresponding 
mean trajectories are also shown (average over 5000 trajectories). First we note the high efficiency, for the T = 0.25 
case. About 95% of the hydrolysis events lead to forward motion. Such a high efficiency is seen in real motor proteins 
at least at low ATP concentrations [tL4|Jl5f| . 

Coming back to the drag force we note that its form has to fulfill certain conditions. Firstly it should act through 
pushing the system into a particular "corner" of configuration space. Secondly, the drag force should not change the 
ATP parameter A directly. However, the reaction rate, set by the term — [ATPj^o, can be influenced through forcing 
the dynamical variables tpa, (p% and P2 ■ Thirdly, the drag force influences the motor function through changing the 
rate of entering the state (<pi, ^2) = (1, 1) — * (0, 0) vs. the (ipi, ip%) = (0, 0) — » (1, 1) rate. An explicit realization that 
fulfills these conditions are 

F(ip , ip u p 2 ) = +f(<po + PaPi + <Pq<Pi<Pz) ■ (9) 

This is not the only possible choice, but given the structure of the Hamiltonian (^), this is a natural one. Another 
possible choice is simply F = f(<po + p\ + ip%). Note that with / positive, the drag force acts against the direction of 
the motion of the motor. 

In figure || we show the efficiency p versus drag / for different temperatures T. We observe that p decreases as the 
drag force is increased. For high temperature the efficiency p decreases approximately approximately linearly with 
/. At lower temperatures one observes that a very high efficiency is maintained even for moderate drag. At larger 
drag, the motor stops moving forward. The stalling at / = 1 occurs because the energy-minimum of Hq, Eq. (Q), 
which is equal to -3, is destroyed by the contribution of the drag term, which is 3/ when / = 1. For / > 1, the state 
(tpo, (pi, (P2) = (0, 0, 0) has lower energy that the state (1, 1, 1), and the ATP-cycling looses its effect. The net result 
is that the motor stalls. The question now is, what temperature should we choose in the model to compare with 
actual motors. The ATP cycle is known to generate about 25kT room and the total binding energy af a protein domain 
is of the same order. Further motor proteins stall at a drag corresponding to a work per step of about !QkT room . 
In the model the total binding energy is 3 and the ATP cycle generates 4 units of energy. Thus room temperature 
corresponds to a temperature in the model of about 0.2. In the discussion above, this is what we mean with low 
temperature. Thus our model predicts high efficiency maintained nearly up to stalling drag. 
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Finally we would like to remark that if we pull on the motor, i.e., employ a negative /, the success rate increases 
for some range of pulls, until finally the efficiency of the motor drops dramatically. This drop is because for / large 
and negative both values of A lead to the same ground state. To be quantitative for the case A = —4: For [ATP] = 
the ground state is (1, 1, 1) for all / < 0. For [ATP] = —4, then for / > — 1/3 the ground state is (0, 0, 0) whereas for 
/ < —1/3 the ground state is again (1, 1, 1). Thus for /> —1/3 the motor switch between the ground states (0, 0, 0) 
and (1, 1, 1) when the ATP variable cycles. However when / < —1/3, the (1, 1, 1) state becomes the ground state for 
both [ATP] values, and the ATP-cycling no longer forces the motor to switch between states. The result is a drop in 
efficiency. In Figure ^[ the drop occurs at a smaller / than —1/3, but seems to approach this value as the temperature 
is decreased. 

We now discuss how the present model is different from existing ones. It has previously been proposed that motors 
could work as ratchets |^-^]. Our model is indeed a type of ratchet, it cycles between two states, one of random 
diffusion and one of directed motion, much like the sawtooth ratchet originally proposed by Prost et al. || . However 
a ratchet mechanism similar to the one studied by Prost et al. would, dependent on the asymmetry of the potential, 
at most make one positive move at half the ATP hydrolysis events. The other half of the ATP cycles, it will remain 
in same potential well, or maybe even stop. Thus the maximum efficiency would be 1/2. In this letter we instead 
propose, that forward motion depends on the relative ordering of two events associated to two different degrees of 
freedom. This results in a motor which is close to 100% efficient. 

In summa the present paper discusses how ordering of events may be utilized to construct a simple motor. The 
motor contains three variables (po, ip\ and ip2- However, the hierarchical structure of the Hamiltonian (Jfy is easily 
extended to any number of variables — and thereby levels. The hierarchical nature of the Hamiltonian ensures that 
each level controls all subsequent levels, that is we have an explicit realization of local control of global motion as 
seen in allosteric proteins. Furthermore, the model suggests a connection between two seemingly unrelated properties: 
folding and motor action. 

FIG. 1. Position of motor vs. time (ATP clock) for two different temperatures: T = 0.25 (upper two curves) and T = 5 
(lower two curves). In each case, both averages over 5000 samples and a single realiztion are shown. The value of the ATP 
parameter A was equal to -4. There was no external drag force acting on the motor. 



FIG. 2. Speed of motor (measured in terms of ATP cycling) vs. drag force / averaged over 5000 samples for different 
temperatures. The ATP parameter A is set to -4. 
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